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We investigate escape processes from metastable states that are driven by non-Gaussian noise. 
Using a path integral approach, we define a weak-noise scaling limit that identifies optimal escape 
paths as minima of a stochastic action, while retaining the infinite hierarchy of noise cumulants. This 
enables us to investigate the effect of different noise amplitude distributions. We find generically a 
reduced effective potential barrier but also fundamental differences, particularly for the limit when 
the non-Gaussian noise pulses are relatively slow. Here we identify a class of amplitude distributions 
that can induce a single-jump escape from the potential well. Our results highlight that higher-order 
noise cumulants crucially influence escape behaviour even in the weak-noise limit. 


Escape from a metastable state underlies a large vari¬ 
ety of phenomena in physics, chemistry, and biology [T] 
and has accordingly received an extraordinary amount 
of attention since Kramers’ seminal work [2j. Recently, 
the problem has seen renewed interest in continuous sys¬ 
tems driven by non-Gaussian noise. Poissonian shot noise 
(PSN) has been investigated in the context of full count¬ 
ing statistics EMU, where the escape rate in e.g. Joseph- 
son junctions can be related to the current statistics. In 
the limit of weak noise intensity D , the Kramers rate 
assumes the exponential form r = e~ s,/ D , generalizing 
the classical Arrhenius result. The key challenge is to 
determine the action S in the presence of the infinite hi¬ 
erarchy of noise cumulants that describes non-Gaussian 
noise. Previous studies have truncated this hierarchy, 
typically after the cubic term [3® HP- Analytical results 
retaining the full hierarchy are available only for particu¬ 
lar amplitude distributions, such as one-sided PSN with 
exponentially distributed HM3 or constant amplitudes 
EOS CEB]. On the other hand, the exponential scaling 
in 1/D breaks down for escape driven by Levy-stable 
noise (LSN), as has been established both in simulations 
[MMI and analytical work f27l - f?i7j . Here, we present a 
unified framework where exact results for the Kramers 
rate are obtained for a large class of non-Gaussian noise 
types. Our approach characterizes in particular the dy¬ 
namics of the escape process and explains the difference 
in the scalings of PSN and LSN driven escape. 

We consider the time evolution of a single degree of 
freedom q under the influence of a conservative force, 
with potential V, and noise £, q(t) = —V'(q) + £(f). We 
specify the noise £ very generally as the derivative of 
a process with stationary increments. It can then be 
described via the characteristic functional 

Gi=[g] = (e^o ds C(s)s(s)) = e f/<is^( g (s))^ (i) 

where the characteristic function if> has the form [3Dj 

4’(g) = iga - Dg 2 / 2 + A <j>{g) (2) 

(/(g) = J dA ( elsA ~ l 9 A ~ x ) p( a )- ( 3 ) 


The first two terms in Eq. ([2]) represent a constant drift 
a and Gaussian noise of intensity D , capturing the ef¬ 
fect of continuous fluctuations. The term </> is the con¬ 
tribution of discontinuous jumps. This becomes evident 
when p is interpreted as the density of independently 
and identically distributed jump amplitudes Aj and A as 
the rate parameter of a Poisson process generating jump 
times tj. G^ is then just the characteristic function of 
Gaussian white noise with drift and superimposed zero- 
mean PSN z(s) — {z(s)), where z(s) = AjS(s — tj) 
and n is Poisson distributed with mean At. If p(A) is 
not normalizable, e.g. because of a power law divergence 
p(A) oc |A|~“ _1 for small A [31] with 0 < a < 2, it can¬ 
not strictly be interpreted as an amplitude distribution. 
Nevertheless, A p(A)dA still gives the rate of jumps with 
size in the range [A, A + dA]. 

For our weak-noise limit we need to have a scale for 
noise amplitudes, which we achieve by writing p(A) = 
po{A/A q)/Aq in terms of a base distribution po(x). For 
<j) this scaling translates into <f>(g) = <po(Aog), where <f>o 
is given by Eq. © with p replaced by p 0 . Table [i] lists 
the amplitude distributions po(x) we consider, normal¬ 
ized so that ( x 2 ) = / dx po(x)x 2 = 1. We also include 
the case where the cumulants of the non-Gaussian noise 
are truncated after some low order by expanding <j>{g). 
In the following we focus on symmetric noise with zero 
drift (a = 0). A straightforward calculation from G^ 
then shows that the noise covariance is Cov(£(f), £(t')) = 
(D + — t'), where Dj, = A (A 2 ) = AAq. 

The starting point for our approach is a path integral 
for the propagator, i.e. the probability of arriving at qb at 
time t when starting from q a . We adopt an Ito conven¬ 
tion and express the propagator as an inverse functional 
Fourier transform, following [32] : 


P(qb,t\q a ) = 


r(qb,t) 


'(9a,0) 


V[q\ / V 


g_ 

.27r. 


0 -S[q,g] 


( 4 ) 


The action functional is given by S[q,g] = 
f 0 ds £ (q(s), g(s)) with the underlying Lagrangian 
£(q, g) = ig (q + V'(q)) + Dg 2 /2 — A <j>(g). Now introduce 
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Distribution 

Po Or) 

cf>o(k) = 4>a{—ik) = (e kx — kx — 1) 

Gaussian 

e~ x ^ 2 /V2n 

(j>o,G{k) = e h / 2 — 1 

Constant modulus 

[5(* — 1) + S(x + l)]/2 

0o, const (fe) = cosh(fc) - 1 

Gamma (a < 2) 

|a;| — 1— “e — l a3 l/[2r(2 — a)] 

0o , a (fc) = [(1 - k) a + (1 + k) a - 2]/[2 a(a - 1)] 

Exponential (a = —1) 

e“ w /4 

0o, exp(fc) = fe 2 /[2(l - fc 2 )] 

Truncated (j> 


0O,trunc(^) = 2 ^ 


TABLE I: Symmetric noise amplitude base distributions po studied. For Gamma noise, po{x) can be normalized to f dx po(x) = 
1 for a < 0, with the exponential a special case, but not for a > 0. At generic amplitude scale Ao one has 0(fc ) = (po(Aok). 


the scaling g —> g/D , so that C = £/ D with 

£(<7,g) = *9 (<2 + ^'(<?)) +<? 2 /2 - \D(j){g/D). (5) 

In the absence of the </>-component (A = 0), C is indepen¬ 
dent of D and we have the familiar scaling form of the 
Gaussian action. This means that taking D —> 0 yields 
the usual weak-noise limit of the path integral. For A ^ 0 
we need to apply a suitable scaling for 0 . Expanding the 
exponential in Eq. © and using 0(g) = 0o(gA o ) yields 

xfl\ = A o( x2 ) (^g ) 2 ^(s 3 ) (Ag ) 3 

® \D J D 2 2! D 3 3! W 


We now consider the scaling A -A A /D^, A a —> 
AqD", which makes the 0{g n ) term of XD<p scale as 
£)i-/i+n(i/-i)_ exponents /x, v define different scal¬ 
ing regimes as D —> 0, as shown in the inset of Fig. |T| 
In regime I, all orders (n > 2) in g diverge as D —> 0. 
In regime II, there are always some higher orders that 
diverge as D —► 0, while in regime III all orders scale to 
zero as D —► 0 so that one effectively recovers the case 
A = 0. For the particular combination v = \ {ji + 1) with 
/x > 1 (red line) only the g 2 term remains in Eq. |5| as 
D —> 0. The non-Gaussian noise strength D$ oc D —> 0 
here, so this is a valid weak noise-limit but one that re¬ 
duces to effective Gaussian noise. Only for /x = v = 1 
do all orders in g remain in Eq. <[5j) as D —> 0. This is 
the scaling we adopt: it represents a genuine weak-noise 
limit of our generic noise, since D$ oc D —> 0 while the 
infinite hierarchy of noise cumulants is retained. 

After scaling, we obtain for the path integral 


p(qb,t\q a ) 



V 


9 

2t tD 


e -S[q,g]/D , 


with Lagrangian 

£(g,g) = *g(g + v'{q)) +g 2 /2 - (8) 

and <j>(g) = ^o(g^o)- We now drop the tildes as only 
rescaled variables are used below. As D A 0 we can 
calculate p{qb,t\q a ) within a saddle-point approach as 
p(qb 1 t\9a) — where = indicates the asymp¬ 

totic behaviour for small D. The optimal paths q *, g* 
are determined from the condition SS = 0, which leads 
to the coupled Euler-Lagrange (EL) equations 


q = -V'(q)+ig-i\(j>'(g) (9) 

9 = V'\q)g, (10) 



FIG. 1: (Colors online) Optimal excitation paths for the dif¬ 
ferent 4>o (see Table[Ij) in the potential V ( x ) = x 4 —6x 2 —2x+5, 
where q a = —1.64, qb = —0.17. Parameter values: A = 1, 
Ao = 1, a = 1.2, 6=1/2 (corresponding to fourth order trun¬ 
cation of </>ex P ). Thin lines: paths directly sampled from the 
weight Q with g integrated out, for D = 0.01. (a) Mean sam¬ 
ple path compared with theory for <f> a . (b) Different scaling 
regimes of non-Gaussian noise given by the scaling exponents 
M> v - 


The boundary conditions are q( 0) = q ai q(t) = qb, and 
accordingly we write the optimal action as S(qb,t\q a )- 
Comparing Eq. © with Eq. ([8]) shows that it can be ex¬ 
pressed as a functional of g* ( t ) only. On physical grounds 
we require real solutions for q* and hence purely imagi¬ 
nary g*, so we define k(s) = ig*(s) and get for the action 


S(qb,t',q a )= / ds 
Jo 


) k 2 + A ( k(j>'(k ) — 4>{k)) 


( 11 ) 


where </>(fc) = <f>(g) = <j>(—ik). Evaluating <p(g) for imagi¬ 
nary arguments in this way requires a continuation into 
the complex plane and so a nonzero radius of convergence 
for <f>(g), a condition that is met if p(A) has tails that are 
no heavier than exponential. One could formulate the 
solution of ( M ) in terms of a Hamiltonian U with the 
conjugate momentum dC/dq = ig: C = igq — T-L(q,ig). 
The EL-equations are then just Hamilton’s equations de¬ 
scribing motion with a conserved Ti. 

Let us now consider specifically escape from a 
metastable state q a , located at the minimum of the 
metastable basin of V, to the top of the potential barrier 
at qb (> g a )- For Gaussian noise, the path integral solu- 
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tion of this problem is essentially analogous to the quan¬ 
tum mechanical tunneling problem treated in a semiclas- 
sical approximation [33] (see also [341136] for a discussion 
in the statistical mechanics context). The theory of large 
deviations gives the dominant scaling of the escape rate 
r in the weak-noise limit as EH [35] 

r 2* e -S{q b \q a )/D ^ S(q b ;q a ) = inf S(q b ,t;q a ). (12) 


The optimal path that provides the lowest (infinimum in 
Eq. @) action is achieved for t —> oo [33]. For A = 
0, S(q b ;q a ) can be determined analytically as twice the 
height of the energy barrier, AV = V(q b ) — V(q a )- Since 
D = 2 T for thermal noise, Eq. (12) thus recovers the 
Arrhenius result [lj. In this case, the optimal escape or 
“excitation” path is the time-reverse of a deterministic 
relaxation path from q b to q a [31H3S] . 

For A ^ 0, deterministic relaxations with g* = 0 still 
solve the EL equations and have zero action, but their 
time-reversal no longer gives the excitation paths. This 
is clear from the predictions for different amplitude dis¬ 
tributions in Fig. [TJ which we have confirmed by direct 
path sampling. The optimal escape paths have the char¬ 
acteristic instanton shape: for large t the system spends 
most of its time close to q a and q b while the actual bar¬ 
rier transition is sharply localized in time. The key ob¬ 
servation is that the instanton shape varies with </>, while 
the deterministic relaxation and its time reverse are en¬ 
tirely independent of <f> and A. Moreover, we observe 
that the optimal action S(q b ;q a ) is reduced compared 
to the Gaussian limit of the noise for a range of small 
A values: the non-Gaussian noise makes escape faster. 
Differences between amplitude distributions become pro¬ 
nounced especially in the limit A —> 0: the Gaussian case 
is approached continuously with the low-order truncated 
<fi and for constant and Gaussian distributed noise am¬ 
plitudes, though the approach is extremely slow for the 
latter (Fig. [2]). For exponential and Gamma noise, the 
action is discontinuous at A = 0: as A —0 it converges 
to a value considerably smaller than 2AV. Puzzlingly, 
for a > 0 the small A regime appears inaccessible, with 
q* becoming complex below some threshold. 

In order to understand these surprising observations, 
we proceed analytically and integrate out g directly from 
Eq-0. In the weak noise-limit this can again be done by 
saddle point integration and gives an action for q alone. 
(Technically we discretize into small time intervals dt and 
take dt —> 0 after D 0.) Defining /3(fc) = k 2 /2 + Xcj)(k), 
the resulting action is S[q] = f dsn(q(s) + V'(q(s))). 
Here 7r(-) is the Legendre transform of /?(•), i.e. n(f) = 
max[fc/ — /3(k)] with the maximum taken over the range 
of k where <p(k) remains non-singular. Note that the 
function 7r is not equivalent to the Hamiltonian H since 
we have integrated out the momenta ig. In fact, since 
q+V 1 (q) = £ from the original equation of motion, the ac¬ 
tion S'jg] = f ds 7r(£(s)) gives the weight of any trajectory 



In A 

FIG. 2: (Colors online) The normalized action Snorm = 
S(q b ; q a )/ (2AF/(1 + AAjj)) with Ao = 1 obtained for the 
same potential as in Fig. [l] The reference value 2AV/(1 + 
AAq) is the action for the Gaussian limit of <j>, i.e. truncated 
after the quadratic term. The (j> a curve is dashed where op¬ 
timal paths cannot be found from the EL equations because 
they involve jumps. Inset: S n orm as a function of InAo for 
0ex P and 4> a , at fixed = AAj) = 1. For large Ao (small A) 
the curves converge to So = / dq min(2F'(q), 1/A 0 ). 


of the noise (averaged again over small dt) in the large 
deviation limit D — > 0. The function 7 t(£) thus general¬ 
izes the simple quadratic £ 2 /2 appearing in the Wiener 
measure exp[— J ds £(s) 2 /(2D)] for Gaussian noise. 

One can now think of q(t) as a path in the ( q , 7))-plane, 
with v = q. Then the action reads S = f dqn(v + 
V'(q))/\v\ and for each q we can find v = q simply as the 
minimum of ir(v + V'(q))/\v\. We do not need to enforce 
the total time constraint t = j dq/\v\ as we want t —> oo 
and the integral automatically diverges at both ends for 
paths between stationary points of V. The trivial global 
minimum is v = —V'{q). which describes deterministic 
relaxation. For an excitation from q a to q b > q a we need 
77 > 0, on the other hand. The condition for a minimum 
of 7 t(t7 + V'(q))/v with respect to v for excitation paths 
can be cast in the form 


V\q) = m/k, 

v = p{k) - V'{q) 


(13) 

(14) 


using basic properties of Legendre transforms. Here, k 


has to be found from Eq. (13) and then gives v using 


Eq. ( |14[ . This implicitly defines a function v = q = 
E(V'{q)) and hence characterizes the shape of the exci¬ 
tation path. Moreover, we obtain the action simply as 


S = 


d qk{q), 


(15) 


Eqs. (13-15) re- 


where k(q) is the solution of Eq. 
produce existing results for special cases. In the Gaus¬ 
sian case (A = 0) one has (3(k) = k 2 / 2; thus k(q) = 
2 V'(q) and v = V'(q ), the expected time reverse of 
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the relaxation path. For escape driven by one-sided 
PSN without a Gaussian component, as investigated in 
Refs. [ 12 l [T 8 l fT 6 l fT 8 ] . we have j3(k) = XA^k 2 /[2(1 —A 0 k)\ 
for exponential and /3(h) = X (e A ° k — A 0 k — l) for con¬ 
stant amplitudes. This yields k(q) = 2V'(q)/(X^ j + 
2 A 0 V'(q)) H2»II31IIS] for the exponential case and k(q) 
as the solution of k = ln(l + k(A 0 + V'(q)/X))/A a [TS] 
for the constant case, respectively. 


To use Eq. (13) to understand the observations in 
Fig. [ 2 ] we distinguish three cases: (a) </> has no singular¬ 
ities for real k, which includes our scenarios </>truno 4>G 1 
^const; (b) 4> has singularities for real k beyond which 
it is not defined, and diverges as these singularities are 
approached; this happens for </> exp or more generally <f) a 
with a < 0, with singularities at k± = ±1/A 0 ; (c) (j> has 
singularities but remains bounded on approaching them, 
as for our (j) a with a > 0. Now write Eq. (13) as 


A = 


(j)(k) 


(V'(q) - k/2) = R(k ) 


(16) 


The RHS R(k) diverges as 1/k for k —> 0 and decreases 
as k increases. In case (a) defined above, it then hits 
zero at k = kc = 2V’(q). The limit A —>■ 0 thus gives 
back purely Gaussian behaviour (A = 0) as we found. 
For case (b) the same is true if V'(q) is small enough 
for kc to lie in the allowed range [k-,k+]. For larger 
V'(q), R(k ) goes to zero already as k —>■ k + because </>(&) 
diverges, so that in the limit A —> 0 we get the non- 
Gaussian solution k = k + = 1/Aq. As then also <f>'(k) 
diverges, the corresponding trajectory q*(s ) must have 
a section with infinite slope, i.e. a jump. For A —> 0, 
k(q) has the value kc = 2 V'(q) or k + = 1 /A 0 depend¬ 


ing on which is smaller, giving for the action (15): S 


S 0 = f dgmin( 2 E'(< 7 ), 1/A 0 ). This is clearly lower than 
the Gaussian limit S = 2AVA Therefore, non-Gaussian 
noise of infinitesimal rate A discontinuously lowers the 
effective barrier for escape. This effect occurs for large 
enough noise amplitudes, specifically 2A 0 y^ ax > 1 where 
V^ ax = maxF'(g) for q £ [q a ,Qb\- We emphasize that 
truncation of the cumulant hierarchy at any finite order, 
which gives polynomial 0, always brings us back to case 
(a) and cannot reproduce even qualitatively the discon¬ 
tinuous behaviour for A —> 0. 

In case (c), small A requires k —► ka due to the bound¬ 
edness of cf>. However, if k+ < kc, then this point is 
outside the allowed range of k , and so for small enough 


A there is no solution for k in Eq. (13), in agreement 


with our observation for Gamma noise with a > 0 in 
Fig. jij Small A here means A < (V'(q) — Ar_|_/2), or 

equivalently V'(q) > Vq = (3(k + )/k + . In the range of q 
where this is true, one can now check that tt(v + V'(q))/v 
is monotonically decreasing for v > 0 , reaching the 
limit k. |_ for v —> 00 : the optimal velocity is infinite, 
^(^'(g)) = 00 , and there must be a jump in the optimal 
path. To the action this then contributes k + f d q where 


the integral covers the g-range of the jump. Our con¬ 
clusion is therefore that for case (c) the optimal paths 
will have jumps even for nonzero noise rates A, provided 
that V^ ax > kg'; for Gamma noise with a > 0 one has 
Vq = {1 + AAq(2“ — 2)/[a(a — 1)]}/(2A 0 ), which as in case 
(b) scales with 1 /A 0 if we keep D$ = AAq constant. The 
presence of a jump in certain parameter regimes indicates 
that there is competition between escape via a sequence 
of small amplitude noise steps, and waiting for an atyp- 
ically large ( 0 ( 1 )) noise pulse to kick the system close 
to the barrier. The latter is preferred when the noise 
amplitudes have an exponential tail with large enough 
scale A 0 , and when non-Gaussian noise pulses arrive at 
a sufficiently small rate to make successful accumulation 
of many small pulses less likely. 

One can push the analysis further, e.g. to show that 
in the limit A — > 0 at constant Aq the function S(-) and 
hence the excitation paths become identical for Gamma 
and exponential amplitude distributions, i.e. independent 
of the exponent a. The limiting shape is S(E'(g)) = 
V'(q) for V’(q) < 1/(2A 0 ), and S(V'(q)) = 00 other¬ 
wise; thus the optimal excitation paths consist of ini¬ 
tial and final segments of time-reversed relaxations, con¬ 
nected by a jump, and the resulting action is just So- 
For large A 0 one then has So ~ (qb — q a )/Ao and the 
excitation path becomes a jump directly from q a to qb- 
In the limit A 0 — 00 our Gamma noise retrieves LSN 
where p(A) oc |A|” 1_Q . Our results show that the effec¬ 
tive barrier governing the escape rate then vanishes, and 
hence clarify why the exponential scaling with 1/D must 
break down in this case as found both numerically |19l - 
[2B] and analytically [27H25] . Note that our arguments for 


all cases (a,b,c) above imply that the solution of Eq. (16) 


always obeys k < kc = 2 V 7 (< 7 ), and hence the action (15) 
is lowered by the non-Gaussian noise, S < 2/S.V. Fig. |2] 
suggests that even S < 2AE/(1 + AAq), i.e. the effective 
barrier is smaller than for the equivalent Gaussian noise. 

In summary, we have defined a weak noise scaling limit 
that allows the effects of non-Gaussian noise to be fully 
accounted for. We have identified a class of amplitude 
distributions that lead to optimal escape paths contain¬ 
ing jumps; these encompass the entire path for large 
noise amplitudes. We have shown that even very rare 
non-Gaussian noise (A —> 0) can discontinuously lower 
the effective barrier for escape, an effect that cannot be 
captured by approximations that truncate the hierarchy 
of noise cumulants. It would be highly interesting to 
observe such non-Gaussian escape events in an experi¬ 
mental system, e.g. in Josephson junctions [551 HD] . We 
remark that our path integral approach can be extended 


to give the prefactor in Eq. (12), by including the fluctu¬ 


ation determinant and regularizing multi-instanton con¬ 
tributions [34, [35j [4T|. Application of our approach to 
systems with more than one degree of freedom could also 
answer open questions such as the effect of non-Gaussian 
noise on the selection of transition configurations [42j . 
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